Exact Solution of Kemmer Equation for Coulomb Potential 
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Abstract 



This article illustrates the bound states of Kemmer equation for spin-1 par- 
ticles. The asymptotic, exact and Coulomb field solutions are obtained by 



I. INTRODUCTION 



using action principle. In the conclusion the energy spectrum of spin-1 par- 
ticles moving in a Coulomb potential compared with the energy spectrum of 
spin-0 and spin-1/2 particles. 
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In the present article we analyzed the problem of massive, charged particles of spin-1 in 
a Coulomb potential using action principle. The quantum mechanics of charged, massive, 
spin-1 bosons in the presence of an external field, was studied for many different situations 
using different techniques |H,|||||| ■ These works especially included an investigation of the 
Ph! solution of the equation in the presence of a magnetic field. These techniques are far more 
difficult to employ when the anomalous magnetic moment coupling terms are included into 
the equations of motion. 

Relativistic spin-1 particle was at first described by Kemmer in 1939 ||. Kemmer equa- 
tion is a Dirac type equation, involves matrices obeying a different scheme of commutation 
rules. It is reviwed because of the interest in the quark-antiquark bound state problem. Cor- 
ben and Schwinger were the first to present a general theory for massive particles of spin-1 
in an external field and arbitrary magnetic moment ||. They generalized the equations of 
Proca and studied the motion of such particles in an external Coulomb field. Their work 
provides the basis for many other studies. 

This article is the first that illustrates the exact solution of Kemmer equation for massive 
spin-1 particle in the presence of a Coulomb potential. The method used here is a generally 
accepted method to solve the relativistic particle equations. The Lagrangian density of 
Kemmer Hamiltonian was separated into radial and angular parts by using the spin space 
rotation operators of SO (3) group in the action principle and then the radial equations of 
motion are obtained by using Euler-Lagrange equations of motion. Finally the asymptotic, 
exact and Coulomb field solutions of these equations are discussed. 



*E-mail address: alihavare@mersin.edu.tr 
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II. RADIAL LAGRANGIAN 



The Dirac-like relativistic Kemmer equation for spin-0 and spin-1 particle of mass m is 

(/3%-m)tt(o0 = (1) 
where the 16 x 16 Kemmer matrices f3^ (/i = 0, 1, 2, 3) satisfy the commutation relation 

P»P V P X + (3 X (3 U (3^ = g^ u (3 x + g Xu ^ (2) 

and represented as 

^ = y* ® i + 1 ® y* (3) 

with 7 M usual Dirac matrices. The dynamical state \I/ (a;) is a 16-component and includes 
both spin-0 and spin-1 particles, tt^ is 4-vector electromagnetic momentum. After a basic 
step Eq.([L|) can be written in the form 



(7 ® I + I ® 7 ) tt - ( c? ® 7° + 7° ® c?) 



7717° g) 7 



where a are 4x4 matrices and given in terms of 2 x 2 Pauli matrices as 



(4) 



a: 



7°7> 



a 
7j 



(5) 



The Lagrangian density of this Hamiltonian is 



L (x) = & (7 ® I + I (g> 7 ) 7T - ("a ® 7 + 7 ® ~a ) ■ 7T* - 7/17° <g> 7 ^ (x) . (6) 

This form must be separated into radial and angular parts to obtain the radial equations. 
The approach of separation used here is a technique employed by Barut and Unal ]7| , which 
introduces a rotation aligning the z-axis of the coordinate system with the radial direction 
r. In the identical-two-fermion system spin space the rotation operator is 



S R = e5 



(fia ± 



where a 1 is the first component and cr 3 is the third component of Pauli matrices. 
This operator preserves the wavefunction resulting 



S f 



( -a (r,9,<p) \ 
ib (r, 6, if) 
ic (r, 9, ip) 
V d(r,9,(p) J 



where a,b,c,d are 4-component spinors. Hence, Eq.(^) takes the following form: 

= (-at,-i6t J -ict J dt)s , a 1 [(7 ®I + I®7 )7r - 



(7) 



(8) 



2 



(c? <g> 7° + 7° <g> "a ) • V - rrry <g> 7°]^ 
The 3-space rotation has the following effect in the spin space: 



f-a\ 

ib 



ic 

V d J 



(9) 



Sr 1 



( a e \ 




( al \ 


a* 


S R = 








U 3 ; 



(10) 



Therefore, in the square paranthesis of the Lagrangian a r , ag, are rotated into cui, a 2 , a 3 
respectively. Then Lagrangian density can be rewritten as 



= (-a f , -ib\ -ic\ d f ) {(7 <g> I + I <g> 7 ) tt + $[(7° <8> a 3 + a 3 7°)<9r 
+ i( 7 ° <g> a + (a~ - dJ) + 7° ® (d+ - ^ + ) + « + (V - 9_) <g> 7 



+a~ (d + - tr + ) <g> 7 )] - m7° <g> 7 } 



f-a\ 

ib 
ic 
\ d J 



(11) 



The d + and cL operators in Lagrangian are helicity lowering and raising operators, respec- 
tively and given by 



d + = d e + -r—^d,n + -a 3 cot 9 
sin 6* 2 

6L = -d e + + t^ 3 cot 9 

sin 6* 2 



(12) 



To evaluate the matrix products in Lagrangian we first define new 4x4 quantities in bloch 
form: 



u 1 = a Cg> I 



(13) 



where I is a 2 x 2 unit matrix, and 



a 1 + ia 2 



a = \ (a 1 - 



(14) 



Then we can define 16 x 16 matrices as follows: 
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a ( 2 ) = I <g) a 

7{i)=y i ®I (15) 
7(2) = 1 ® 7 M 

where I is now 4x4 unit matrix. In terms of these we can write the following matrices 

a t) = 2 ( a (i) ±ia (i)) 

«f 2 ) = \ ("(2) ± * a (2)) ( 16 ) 

Straightforward algebra allows us to rewrite the Lagrangian in the form: 

L (V) = {2iXQa ] a - 2n Sd + [aVf 6 + aVfc - ftVfc - tfa%a - b ] a\d 

-c + a\a - c + ald + d + a\b + tf 1 " afc](d r + -) + ^-[a + (~a 1 x 7r 2 ) 3 6- 
a + x (? 2 ) 3 c - 6 1 " ( (?i x "^ 2 ) 3 a + 6 + (Tt'i x d + c + (a^ x ~a<^f a 

—c + (iji x ~a 2 ^) d — d + (ij 1 x ~a^j b + d + x Ir^j c] — - [—a)d + a 2 b— 
—a)d + cr^c + b^d + a 2 a + b^d+a^d + c^d+a^a + c^d + a 2 d — d^d + a^b — d^d + a 2 c] 
+ i[-a t <9_ (7 +6 - atd-atc + tfd-a+a + tfd-afd + Jd-a+a + Jd-a+d 

—d^d-aib — d^d-a^^—ma^a + mtfb + mc^c — md^d} . (17) 

where everything is 16- dimensional. One can verify that the square of total angular momen- 
tum in spherical coordinates 

J 2 = H 2 + cot 6de + + ^} + J 2 , (18) 

can be written in terms of d + and <9_ operators as 

j 2 = hd[d + + dld_] + j 2 r (19) 

This allows us to write the 4-component a,b,c,d bispinors in terms of angular functions: 



«=~£(2J + 1) 

j,m 



( a+ (j, m; r) D+i tm (0, ip) \ 
a (j, m; r) £^ j7n (0, y>) 
a 5 (j, m; r) Z^ m (0, y>) 
V a_ (j, m; r) Di x (5, y>) / 



(20) 



],m 



1) 



/ &+ (j, m; r) £>^ 1)m (0, (p) \ 
b (j, m; r) D J 0m (6, ip) 
h (j, m; r) D 3 d m (6, ip) 
\b- (j, m- r) Di lm (9, <p) J 



(21) 



An 



E (2J + 1) 



3,m 



( c+ {j, m; r) D J +l m {6, (p) \ 
c (j, m; r) D 3 Q m (6, ip) 
Co (j, m; r) D\ m (6, ip) 
(j,m;r) Dl lm (9,<p) J 



(22) 



3,m 



( d + (j, m; r) D J +1>m (6, p) \ 
d (j, m; r) D° %m (6, ip) 
d 5 (j, m; r) D ] d m (6, ip) 
\ gL (j, m; r) DL l m (6, ip) J 



(23) 



where D 3 Xm functions are the transformation matrix elements of a general rotation. A 
+1, 0, 5, — 1 indices of D operator are indicating the spin orientations of particle system. 
The orthogonality relation between these operators is given as [H 



■2jt 



sin 6d6D(l (e, ( p)D{] m (e, ( p)dedip 



2 •' _|_ °mm 



and the effect of helicity raising and helicity lowering operators on D operator: 



d.D J 



d + D 3 



X,m 
j 



j (j + 1)D> 



Ij (j + l)D{_ hm 

Next by calculating both sides separately in Lagrangian it is seen that 



and 



J (J + 1)D a k 



J U + l ) D ^ 



for k — 1,2. 

For Lagrangian density the action is given in the form: 



A 



dt r 2 dr E * f (r, t) D% m . (r) [i^ - H] D{] m (f) * (r, t) 

jj'mm' 



(24) 



(25) 



(26) 



(27) 



(28) 
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where 



dr = sin 



(29) 



We can carry out all the angular integrations with the aid of Eq.(^4j) , thereby obtaining an 
effective radial Lagrangian: 



L (V) = 4vr 2 J] (2j + 1) r [a f , 6 t , c + , d*] 



3,m 



'o - rn) 


T 2 


Ti 





-T 2 


rn 





-T x 


-T x 







-T 2 





T-v 




-(2tt + 



xr 



c 



(30) 



where 



T x = afd r 



T 2 = - 



— O" 1 X °" 2 

2r 



ff 2 " + ^l ff l X °" 2 



(31) 



(32) 



By using the internal transformation matrices which are given as follows 



7T 



v/2 



^ 1 1^ 
110 
0-110 

V -l o 01/ 



7T 



v/2 



/" 1 o -i\ 

1-10 
0+110 

V+i o o i / 



(33) 



we define new quantities in terms of linear combinations of the original spinor components: 

(34) 



7T 



fa \ 






f a + d\ 






\ 


b 


r = 


r 


b + c 








c 


72 


c — b 




^ c 




\d J 






\d — a J 






) 



r ( at fet c t rft ) = _!_ ( ( a t + rft) (fct + c t) ( c t _ 6 t) 

= ( *1 *t ) 

Since, 7T7rt = 7T% =1 we can rewrite the Lagrangian in the form 



dt - a t 



)) 



(35) 
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L (V) = 4tt 2 (2j + 1) r [a f , b\ c f , d f ] ttV 





/ (27T - m) 




Ti 









fa\ 


7T^7T 


-T 2 







-2i 




x ix^ixr 


b 


-Ti 





m 


-T 2 




c 




V o 




T 2 


-(2tt + 


m) / 







(36) 



Then, with the aid of identities given in Eqs.(|4]), (|35D we obtain the radial Lagrangian in 
the form: 



L (V) = 4tt 2 (2j + 1) + (Ti + T 2 ) ^ B - 2n ^ D ] + tf B [m* B - (T, - T 2 ) * A ] 

+¥ c [m^ c + (Ti - T 2 ) tf D ] + *J,[-2tt *a + 0*2 - T x ) V c ~ m^f D ]} (37) 



III. RADIAL EQUATIONS 

The Euler-Lagrange equations of motion are 





dL_ 



(38) 



where q = V A , * B> * c , V D . 

These yield four sets of 4-component equations, ultimately a set of sixteen equations: 



m^ A - (Ti + T 2 ) ^ B + 2n ^ D = 
m^B- (T 1 + T 2 )^ A = 
mm c + (T x - T 2 ) tf D = 
+ (T x - T 2 ) ^ c + 2vr ^A = 



The explicit forms of (T\ + T 2 ) and (Ti — T 2 ) are given as 



(2d r -f \ 



Ti+T 2 



0-- 



V o 



A 

^00- 

A A 



r r 



A 

25 r y 



T x -T 2 



/ £ 

2a 



-7 X 

2 A 



A _2 _ 2(9 A 

r r ' r 

V o £ o / 



where A = (j + 1). 

Then, by defining ^l/^, ^c-, spinors as 



(^ai\ 












(*Dl\ 


^A2 
^A3 


, ^B = 


^B2 
^B3 




^C2 




*D3 

















(39) 
(40) 
(41) 
(42) 



(43) 



we obtain sixteen radial equations from Eqs. (p9|) , ([42|): 
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r r 

m^ A2 - -V B1 + -^ B4 + 2ir ^ D2 = 
r r 

A A 

m^ A3 - -* B1 + -^ B4 + 27To^ D3 = 
m^ A4 - -*B2 - -^B3 + 29 r ^ B4 + 2vr ^D4 = 

A A 

m^Bl - 2d r ^Al + ~^A2 + ~^A3 = 

r r 

m ^ B2 + _^ A4 = o 

A A 

m^B3 Vai + ~^A4 = 

r r 

m^ Bi - -^ A2 - -^ A3 + 2<9 r ^A4 = 
m^ cl + ^ m - ^ D3 = 
m^ C 2 - + 2<9 r ^ D2 + -^ D3 - -^ D4 = 

ry f f 

m^C3 + -*r»i - 2<9 r ^ D3 - -^ D2 + -^D4 = 

ry ry ry 

A A 
m^ C4 + -m D2 ^ D3 = 

r r 

mty D1 + -^C2 ~ ~^C3 + 2vr ^Ai = 
r r 

$ C1 + 2<9 r ^ C 2 + -^C3 + 27T ^A2 = 

r r r 

A „ 2 A 

m^ D3 + -*ci - 2<9 r ^ C 3 - -^C2 + -^C4 + 27T ^A3 = 

ry ry 

m^ m + - ^* C3 + 2vr ^A4 = 



(44) 
(45) 
(46) 
(47) 
(48) 
(49) 
(50) 
(51) 
(52) 
(53) 
(54) 
(55) 
(56) 
(57) 
(58) 
(59) 



In these equations to eliminate the equations that correspond to spin-0 we must analysis 
the spin orientations of spinor wavef unctions. The spinor wavefunction of the system must 
be written in terms of upper and lower components for particle system: 



# = 



Si 



L 2 

s 2 



( L X L 2 \ 
L\S 2 
S\L 2 



where a, b, c, d spinors are given in the form 



L2, 
L 2 _ 



(a\ 
b 

c 



V S 1 S 2 / 


\ 


d) 


(L l+ L 2+ \ 




/ a+ \ 


L\ + L 2 




a 


L\_L 2+ 




a 5 


V 1.1^2 J 




\a-J 



(60) 



(61) 
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d 



s x _ 



S2 + 

s 2 _ 



s 2 _ 



L2, 
L 



2- 



S2 + 

s 2 _ 



/ T Q 

f ivi + D2 + 






(K\ 






— 


bo 


L\_S2 + 






b 5 




/ 




\ o- J 


/ c r 

/ Ol + lJ2 + 


> 






Si + L 2 _ 




— 


Co 








Co 


\S 1 „L 2 _ 






\c-/ 


1 Si + S2 + 


\ 






Sl + S2_ 






do 


Sl_S2 + 






do 


V Si_s 2 _ 


/ 




\d-J 



(62) 



(63) 



(64) 



Since our system is a single-particle system we can omit indices 1 and 2. Hence, from 
Eqs.flfH]) , (0) it is seen that 



and from Eqs. ( |62"D , ( p3| ) 

&o = b b 



a — a 5 



<"■+ 



dn — dr, 



Co = c 5 



(65) 



(66) 



Then, the new forms of ^b, ^c> ^ d spinors are 



•c 





1 a + + d + 








r 


do + do 






^A2 


V2 


do do 






^A3 




\ a_ + d_ 






\*A4/ 




( c+~b+ 


\ 






r 


c - b 






^>C2 




Co - b b 






^C3 




\c--b- 


/ 




V*C4/ 



= — 





'b+ + c+\ 


/ 




r 


b + c 




^B2 


V2 


b 5 + Co 




^B3 




V b- + c_ / 


I 


v*B4/ 




/ d + — a + s 


\ 




r 


do ~ do 








do Q-o 








\ gL — a_ , 


/ 





We can see the following equalities: 

^A2 
^D2 



^A3 
^D3 
^C4 



*S2 
^C2 




*B3 
-^C3 



(69) 



When we use these equalities in Eqs. [44 — 59] we find that some of the solutions are 
identical. By eliminating one of the identitical equations, we obtain 10 equations for spin-1 
particle: 



m^ A i - 2d r Vm + 2-^B2 + 27t ^di = 
r 

m^ A2 tffli + -* B4 + 2vr ^D2 = 



(70) 
(71) 
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m*A4-2-* B2 + 2d r * i? 4 + 27ro*D4 = (72) 
r 

mV m - 2d r ^ A1 + 2-^ A2 = (73) 
r 

m ^ B2 ^r A1 + -m AA = o (74) 

r r 

m* S4 - 2-* A2 + 2d r ^A4 = (75) 
r 

m ^ C2 _|_ 2d r ^ D2 + ~^D2 - ~^D4 = (76) 

r r r 

m^ Dl + 2-^ C 2 + 2TT ^ A1 = (77) 

r 

2 

m^ D2 + 2<9 r ^ C 2 - ~^C2 + 27T *A2 = (78) 

r 

m^ D4 + 2-^ C 2 + 27T ^44 = ( 79) 

r 



IV. ANALYSIS OF THE RADIAL EQUATIONS 

A. Asymptotic Solutions 

If we write the radial equations of spin-1 particle in the limit r — > oo , we obtain the 
following equations: 

m^ A1 - 2d r V m + +2n ^ D1 = (80) 

m^ A2 + 2tc ^d2 = (81) 

m* A4 + 2d r V B4 + 2ir V m = (82) 

m^Bi - 2<9 r ^Ai = (83) 

m^ B2 = (84) 

m^ m + 2d r ^ A4 = (85) 

m^ C 2 + 2d r ^ D2 = (86) 

m& D1 + +27T ^ A i = (87) 

D2 + 2d r ^c2 + 2vro^A2 = (88) 

m# D4 + +2tt ^a4 = (89) 



The solutions of these equations are as follows: 



y A1 = A ie ikr + B ie ~ ikr 


(90) 


m A2 = A 2 e ikr + B 2 e~ ikr 


(91) 


^D2 = A 3 e ikr + B 3 e~ ikr 


(92) 


m A2 = A 4 e ikr + B 4 e~ ikr 


(93) 


* m = A 5 e ikr + B 5 e~ ikr 


(94) 


m m = A 6 e ikr + B 6 e~ ikr 


(95) 
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^C2 


= A 7 e ifcr + B 7 e~ ikr 


(96) 


*B1 


= A 8 e ikr + B s e~ tkr 


(97) 




= A 9 e ikr + B 9 e~ ikr 


(98) 



where 

e = E 2 -— (99) 
4 



B. Exact Solutions 



For free case the exact solutions of spin-1 particle are found in the form 
(p) = l - [ApJ, (p) + Axe* + Bxe"* 
i [-ApJ, (p) + A ie v + S ie 



#A4 (p) 
*B1 (p) 



2 

!{£>[./,■ (p) + ip ( (p) - J, +1 (p))] + ife* + Fe-*} 
(p) = (p) + \p (Jj-i (P) - Jj + i (P))] + ^ P + Fe"^} 



*D1 (P) 
^D4 (P) 



CpJ 3 (p) + Ge*" + He~ ip 



2 



-CpJj (p) + Ge ip + He- 



lp 



^B2 = (P) - ^A4 (P)] 

pm 



where p = kr. 



(100) 
(101) 
(102) 
(103) 
(104) 
(105) 
(106) 



C. Coulomb Field Solutions 



For Coulomb field 



7r = P - e$ = E + 



(107) 



If it is used in spin-1 equations we obtain the following expreses for ^a, ^b, ^c, ^ d spinors: 

1 



^ai (p) 



AR n ,i (p) + CR nX (p) 



^A4 (p) = - [-AR n:l (p) + CR nX (p) 



Vm (p) 



1 / ka 
— IE+ — 
m\ p 



AR ntl (p) + CR nX (p) 



1 / ka\ 
^D4 (p) = - [E + — Ai^ (p) - Ci^ r (p) 
my P / 



(108) 
(109) 
(110) 

(111) 
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where 



*bi (p) = [AR n + (p) + CR nX (p) 
^B4 (p) = ~d p \AR nil (p) - CR nX (p) 
^ B2 = -BR n j (p) 



2 2a T? 

a = e , n — —hi 



(112) 
(113) 
(114) 

(115) 



and 



Rn,i (p) = e ±lp [±2ip\ 



'('+) 1 



^ i^i (^#+l)Ty,l + ^, ±2. P ) 



(116) 



is confluent hypergeometric function. 



V. DISCUSSION 

In relativistic dynamics the exact solutions of the wave equation are very important 
because of the understanding of physics that can be brought by such solutions. These 
solutions are valuable tools in determining the radiative contributions to the energy. 

Exact solutions of the Kemmer equation with potential interaction are very rare. There 
has been a great deal of interest in the Kemmer equation for a constant magnetic field but, 
on the other hand, we first solved exactly the Kemmer equation for relativistic Coulomb 
potential. Since Kemmer equation is a two-body Dirac-like equation it includes both spin-0 
and spin-1 particles. 16 solutions we obtained include both spin-0 and spin-1 particles. The 
identical solutions of these 16 equations represant the spin-0 particle. By eliminating one of 
the identitical equations, we obtained 10 equations for spin-1 particle. 

The asymptotic solutions of spin-1 particle are found in terms of periodic functions. 
These are important in physics in order to leave the components of spinors finite. The free 
solutions of spin-1 particle are found in terms of Bessel's functions. In these solutions the 
(p) , ^ bi (p) , ^di (p) and ^ aa (p) , ^54 (p) , ^D4 (p) can be interpreted as helicity states 
of the particle and anti-particle, respectively. 

For Coulomb potential the solutions are obtained in terms of confluent hypergeometric 
functions: 



R n ,i (P) = [±2.p]H 1+ V <4+) J lFl (1 + + 1)^*1, 1 + yi(l + l), ±2*p) 

where 1(1 + 1) — j (j + 1) — a 2 . 

Here 1F1 (a, b; ±2ip) denotes the degenerate hypergeometric function. For large p this 
functions behaves as 
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(117) 



Demanding that F be normalizable, that is J °° dr |F| 2 = 1 implies that [r (a)] 1 must vanish. 
This is the desired quantization condition M . For our solution the quantization condition is 




-n 



(118) 



which leads to 



E 



± (m/2) 



4 Q 2 



(119) 



Hence, by using the Coulomb field solutions of Kemmer equation we can predict the energy 
levels of spin-1 particles. This spectra is similar to the energy spectrum of spin-| particle 
moving in Coulomb potential which is given by 



E 



±m 



1 + 



(120) 



In a previous work |10j, the bound states of a spinless charged particle in the Coulomb field 
has been obtained. The radial wave functions of spin-0 particle for Coulomb interaction was 
obtained in terms of associated Laguerre polynomials: 



F nJ (p) = Ne-'fiftlff 1 (p) 



(121) 



where n = A — I — 1 , X = n — J + I and N is a normalization constant. The eigen-energy 
of spin-0 particle for Coulomb potential is 



E 



ni 



(122) 



(n-J-iy 



From the energy spectrums of spin-1 and spin-0 particles one can see the first order contri- 
bution (ja 2 ) to the energy that comes from the spin-Coulomb field interaction by expanding 

the + y (j + |) — ct 2 ^ quantum number in series in the case of spin-1. It can also be 

seen how these contributions depend on spin by considering the spin- 1/2 case. 

The exact solutions of spin-1 particle for Coulomb potential is also important in QED. In 
QED the vacuum polarization in Coulomb field is calculated using Coulomb wave functions 
PH . The solutions obtained here can be used for this purpose. 
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